Assessment of density functional approximations: long-range correlations and 

self-interaction effects 
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The complex nature of electron-electron correlations is made manifest in the very simple but 
non-trivial problem of two electrons confined within a sphere. The description of highly non-local 
correlation and self-interaction effects by widely used local and semi-local exchange-correlation en- 
ergy density functionals is shown to be unsatisfactory in most cases. Even the best such functionals 
exhibit significant errors in the Kohn-Sham potentials and density profiles. 
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I. INTRODUCTION 

The Kohn-Sham (KS) Q formulation of Density Func- 
tional Theory (DFT) 0] is in the present day the most 
popular method in electronic structure calculations. In 
this scheme, the exact ground state energy and electron 
density could be found self-consistently if the exchange- 
correlation (XC) energy functional Exc[n] was known. 
Exc{n] contains all the quantum many-body effects of 
the electron system, but very simple mean-field prescrip- 
tions like the Local Density Approximation (LDA) often 
suffice to obtain accurate results for a wide variety of 
systems at an affordable computational cost. 

However, there are several problems that are well be- 
yond the capabilities of the local approximation and any 
semi-local extension thereof, because of the clear mani- 
festation of the very non-local nature of electron-electron 
correlations. For example, the long-ranged van der Waals 
interactions, the image potential at metal surfaces and 
clusters, or several pathological behaviors of the exact 
XC potential cannot be described at all using simple 
XC functional models Nevertheless, these limita- 

tions should not be important in many other situations, 
and new XC functionals are being proposed with the aim 
of reaching chemical accuracy while keeping the imple- 
mentation ease of the local density based approxima- 
tions 0, . This could open the appealing possibility of 
making predictive studies of relevant aspects of Quantum 
Chemistry such as reaction paths, atomization energies, 
and bond lengths and energies. 

The purpose of this paper is to bring further insight 
in the capabilities and limitations of local and semi-local 
XC functionals, showing that even in very simple prob- 
lems the complexity of the quantum electron-electron 
correlations might prevent any of these approaches from 
properly describing the ground state properties of such 
systems. We do not intend to make a comprehensive 
assessment of mcan-ficld-likc approximations, but just 
to provide a representative common picture of all them. 
Hence, among the realm of proposals existing in the liter- 



ature, we have chosen the well-known LDA prescription 
by Perdew and Wang [(|, the Generalized Gradient Ap- 
proximation (GGA) by Perdew-Burke-Erzenhof [?J, and 
the very recent meta-GGA (MGGA) proposed by Tao 
et al Q • These three approaches have the virtue of be- 
ing designed using general considerations (i.e. they do 
not include empirical parameters). Furthermore, they 
can be seen as a coherent set of conceptual progressive 
improvements starting from the strictly local approxi- 
mation, then considering the dependence on the density 
variation, and finally including information from KS or- 
bitals through its associated kinetic energy density. 

We will study a very simple but non-trivial system: two 
electrons confined within a sphere of hard walls, whose 
solution has been found through accurate Configuration 
Interaction calculations 

UGH- This system is an inter- 
esting benchmark reference due to the following reasons. 
First, its simplicity: the density is isotropic, hence hav- 
ing a simple mathematical one-dimensional problem re- 
stricted to the radial coordinate. Second, its ground state 
has singlet spin configuration. As a consequence, Pauli's 
correlation (exchange) between the electrons is absent 
and the exchange energy E% [n] just corrects the spurious 
electron self-interaction in the classical Hartree electro- 
static energy Wh [n]. That is, the Coulomb correlation is 
actually the only source of quantum many-body effects 
in this system. Finally, different correlation regimes can 
be easily achieved by varying the radius R of the confin- 
ing sphere. Thus, at small R (high mean density) we are 
in the low-correlation limit where the confinement by the 
sphere dominates over the electron-electron interaction, 
so having a system with an atomic-like behavior. By in- 
creasing the value of R (decreasing the mean density) we 
gradually enter into a highly correlated regime in which 
the correlation exhibits long-ranged and anisotropic ef- 
fects. 

Then, this simple system offers an excellent scenario 
to assess essential features of functional approximations 
to the XC energy. In particular, since it is impossible 
within the present formulation of semi-local functionals 
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TABLE I: Comparison between the exact exchange and correlation energies for several sphere radii R and the results given by 
the local and semi-local functionals evaluated on the exact density profile. The exact total energy is also included to illustrate 
the increasing importance of correlation for high R. 



R 


Z? ex 

■C'tot 


tti ex 


ttiLDA 


^x 


ttiMGGA 

^x 




ttiLDA 


TTnGGA 

hj C 


ciMGGA 


1 


11.5910 


-1.7581 


-1.5230 


-1.6843 


-1.7805 


-0.0507 


-0.1424 


-0.0775 


-0.0647 


5 


0.7016 


-0.3335 


-0.2914 


-0.3213 


-0.3412 


-0.0383 


-0.0715 


-0.0501 


-0.0421 


10 


0.2381 


-0.1592 


-0.1407 


-0.1550 


-0.1651 


-0.0288 


-0.0481 


-0.0362 


-0.0301 


25 


0.0633 


-0.0590 


-0.0537 


-0.0596 


-0.0634 


-0.0163 


-0.0257 


-0.0201 


-0.0166 


50 


0.0249 


-0.0278 


-0.0262 


-0.0296 


-0.0311 


-0.0093 


-0.0151 


-0.0116 


-0.0095 



to achieve the exact exchange energy for arbitrary one- 
and two- electron systems (— Wh [n] and — Wh [ n ] /2 re- 
spectively), we can easily see how important is this limi- 
tation for two-electron densities with very distinct mean 
densities. On the other hand, LDA and GGA suffer from 
spurious correlation self- interaction, a limitation which is 
corrected by the MGGA Nonetheless, the proper 

self-interaction correction does not guarantee the overall 
accuracy of the correlation functional, and the actual role 
played by non-local correlation effects has to be checked 
carefully. 

As a first test, we will evaluate these functionals over 
the exact densities, i.e. in a non-self-consistent fashion, 
comparing the XC energies and potentials with the ex- 
act ones. Then, we will present the fully self-consistent 
solutions, in such a way that we will assess not only 
the self-consistentent energies but also the DFT densi- 
ties that minimizes the corresponding total energy func- 
tionals. Atomic units (H = m e = e = 1) will be used 
throughout the paper. 



II. RESULTS ON THE EXACT DENSITY 
PROFILE 

The Hamiltonian of the two-electron model system is 
given by 



H 



1 2 



1 



n - r 2 



V(r) 





oo 



i=i 

r < R 
r> R. 



The hard wall described by V(r) impose strict boundary 
conditions at r = R, but does not have a direct contri- 
bution to the total energy functional E [n]. Therefore: 



E[n] = T s [n] + W H [n] + E X c[r 



(1) 



where T$ [n] is the kinetic energy of the fictitious KS non- 
interacting system. For a singlet state, the set of KS 
equations is reduced to a single Schrodinger equation 



^s(r) 



(r) = e0(r) , 



(2) 



where vg = i>h + «x + is the KS effective potential 
which is the sum of the Hartree (j)h), exchange (vx), and 
correlation (vq) potentials. The density is related to the 
ground-state orbital of © through the simple equality 
0(r) = 

Under the exact DFT formulation, if n cx (r) is the 
ground-state density of the system, the corresponding 
eigenvalue e cx must equal the ionization energy E[n ex ] — 
E^\ where i?[n GX ] = E^ t is the two-electron ground 
state energy and E^ = n 2 / (2i? 2 ) is the energy of the 
one-electron system. Thus, the exact correlation poten- 
tial can be written explicitly as 



where we have used the exact relation Ex [n] = 
— Wh [n] /2.|l3 However, it is worth pointing out that © 
is an expression only valid for the exact density profile. 
An exact functional expression for the correlation poten- 
tial vc (r) of an arbitrary spin-unpolarized two-electron 
density n (r) would require to know the external potential 
that defines the two interacting electron system whose 
ground state is n(r) and then include such a potential. 
Then, the simplicity suggested by Q is just apparent. 

The performance of different local and semi-local pre- 
scriptions when evaluating the XC energy on the exact 
density profile is shown in Fig. Q]and Table [I] The LDA 
systematically underestimates the absolute value of the 
exchange energy, whereas the GGA partially corrects this 
trend, although in the low density limit the GGA over- 
estimates | -Ex |- The MGGA behaves reasonably well for 
small radii, which is not a surprise since it reproduces 
exactly the exchange energy of the hydrogen atom and, 
as we said in the Introduction, in this range the model 
system behaves precisely like an atom. Thus, although 
the MGGA does not cancel exactly the spurious Hartree 
self-interactions, it fairly accounts for such a cancella- 
tion in atomic-like systems. When decreasing the mean 
electron density, the system cannot be considered like an 
atom any more and the MGGA greatly overestimates the 
self-interaction corrections to Wh, becoming even worse 
than GGA. 

Regarding the correlation energy Eq [n], the LDA 
shows an evident poor behavior which is improved by the 
GGA although the correlation energies are always much 
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FIG. 1: Percent errors AE = 100 x (_B DFT - E cx ) /£ cx for 
the exchange, correlation, and total energies as functions of 
the sphere radius R. All DFT results are obtained non-self- 
consistently over the exact density. Solid line: LDA; dashed 
line: GGA; dotted line: MGGA; dash-dotted line: EXX+M. 



too negative. On the contrary, the MGGA behaves ex- 
tremely well for all densities. Its relative error in the 
high density limit (around 20%) has a minor influence in 
the total energy, and such an error is less than 5% for 
lower densities. This excellent performance is in agree- 
ment with the conclusions by Seidl et al |l2j| about the 
essentially correct behavior of the meta-GGA correlation 
functional proposed in Ref. 11 (the basic ingredient of 
the MGGA by Tao et al) under uniform scaling to the 
low density limit. We have to bear in mind that in the 
highly correlated regime, the two electrons are localized 
in two different positions. That is, if one of the electrons 
is at a distance r from the center, the probability to find 
the second one is concentrated around a point on the op- 
posite side . Thus, the LDA/GGA main source of error 
in this regime, corrected by the MGGA, is the absence 
of self-interaction corrections. 

In the lower panel of Fig. ^ we present the accuracy 
of the corresponding non-self-consistent DFT total en- 
ergies. The well known compensation of errors between 
exchange and correlation in LDA and GGA is easy to 
observe in the atomic-like limit, but in the high corre- 
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FIG. 2: Exchange and correlation potentials obtained from 
the exact density profile (R = 10). Thick solid line: exact 
results; thin solid line: LDA; dashed line: GGA; dotted line: 
MGGA; dash-dotted line: EXX+M. Note that the shape of 
the LDA XC potential reproduces fairly well the exact one 
excepting in the region around the center of the sphere. This 
error is amplified by the EXX+M prescription, where the ex- 
change part is exactly correct but the correlation potential is 
the same as the MGGA. 

lation regime both functionals fail badly. The MGGA 
does not benefit from any cancellation of errors (it over- 
estimates the absolute value of exchange and correlation 
energies). Hence, it slighty underestimates the total en- 
ergy for small radii, but the error on exchange dominates 
for lower densities where the MGGA perfoms poorly. A 
hyper-GGA designed with the aim of being free of 
any self-interaction error, must yield very accurate total 
energies for all the ranges in our model system. Then, 
its performance should be similar to the presented in the 
lower panel of Fig. ^ where we plot the relative error on 
the total energy evaluated through a hybrid functional, 
EXX+M, where the exchange is calculated exactly and 
the correlation approximated by the MGGA. In this case, 
the error on the total energy is always less than 1 %, and 
for R ~ 1 the absolute deviation is just 7 mHa/e, fairly 
close to the chemical accuracy (around 2 mHa/e). 

From the shape of the potentials it is possible to see 
the underlying physics contained in the approximations. 
The above mentioned localization of the electrons in the 
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highly correlated limit is an obvious consequence of the 
electrostatic repulsion, which tend to dominate over the 
confinement by the wall as we increase the radius of the 
system. Since the Hartree energy contains spurious self- 
interactions, the role of the exchange is to compensate 
partially the effects due to Wn . As we can see in the up- 
per panel of Fig. |3 where vx(f) is plotted for the model 
system with R = 10, the exact exchange potential has a 
mininum in the center of the sphere, favouring an atomic- 
like behavior. The corresponding LDA potential is not 
able to account completely for this exchange attractive 
feature, whereas the overall shift of v x GA is a concomi- 
tant consequence of the local dependence on the density. 
The semi-local potentials exhibit a similar behavior but 
there are unphysical oscillations reflecting the presence 
of the gradient and the Laplacian of the density in the 
expression of the exchange potential. 

On the contrary, the Coulomb correlation enhances the 
localization through a potential barrier located in the 
center of the sphere (see the middle panel of Fig. |5J). 
This barrier reflects a truly non-local correlation effect. 
In fact, for R — 10, the electron density is almost homo- 
geneous around r = 0, and then v}P A (r) is practically 
constant in this region. On the other hand, the LDA 
fits reasonably well the exact potential if r > 5, but this 
partial agreement is completely lost under the GGA and 
MGGA. This overall bad quality of the local and semi- 
local correlation potentials is a general feature in finite 
systems 01 1 although in some cases it could be masked if 
we focus on the total vxc ( r ) ■ As we may see in the lower 
panel of Fig. |2 the XC potentials given by LDA, GGA, 
and MGGA are rather similar (excepting the above men- 
tioned unphysical oscillations). Their shape is close to 
the exact one for r > 5 but, as expected, they do not re- 
produce at all the correlation barrier at r = 0. Under the 
EXX+M prescription, the situation is even worse, since 
the XC potential reaches a minimum at r = due to an 
exact description of exchange which is not compensated 
by an accurate correlation potential. Therefore, although 
the correlation energies given by the MGGA are excel- 
lent, the potentials derived from it are not able to repro- 
duce the non-local effects that manifest themselves in the 
shape of vc( r )- As we will see in the next section, this 
will lead to important deviations from the exact density 
profile when solving self-consistently the KS equation. 



III. SELF CONSISTENT RESULTS 

One of the major advantages of KS-DFT is its fully 
self-consistent character: any previous knowledge of the 
electron density profile is not required excepting for set- 
ting up an initial guess to start the iterative resolution of 
the KS equations. For many purposes, LDA and/or GGA 
give accurate self-consistent densities because the corre- 
sponding approximate potentials are very similar to the 
exact ones in those regions relevant for the calculation of 
the electron density. This justifies the use of more sophis- 
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FIG. 3: Kohn-Sham potential vs (r) for the exact density pro- 
file (R = 10). Thick solid line: exact result; thin solid line: 
LDA; dashed line: GGA; dotted line: MGGA; dash-dotted 
line: EXX+M. In order to allow an easier comparison, the 
potentials have been shifted in such a way that the corre- 
sponding first eigenvalue for each approximate potential is 
equal to zero. The behavior of the potentials near the sphere 
walls is less important due to the role played by the boundary 
condition n(R) — 0. 



ticated functional expressions as a mere correction over 
the self-consistent LDA/GGA densities. Consequently, 
the main effort carried out during the last years had 
been directed towards the improvement of the XC en- 
ergies, paying less attention to the characteristics of the 
XC potential vxc( r )- 

Nonetheless, in the previous section we have seen that 
all the functional considered in this paper fail to re- 
produce the exact XC potential in a region where the 
electron density is far from being negligible. The con- 
sequences can be seen in Fig. [21 where we compare the 
exact KS potential vs(r) with the DFT ones obtained 
from the exact density profile n cx (r). The LDA, GGA, 
and MGGA do not reproduce the exact shape of vg(r) 
and it is reflected by a classical forbidden region greater 
than the actual one: the self-consistent density will be 
pushed towards the walls. The EXX+M model, as co- 
mented, incorporates exactly the attractive character of 
exchange, but the wrong description of Vc(t) makes the 
effective potential less confining than the exact one: this 
hybrid approach tends to concentrate the density around 
r = 0. 

These deviations from the exact two-electron density 
profile, which can be quantified through the expression 



Sn = - I dr |n ox (r) - n(r)| , 



(4) 



should not be important for high mean densities. In this 
atomic-like limit the external confining potential dom- 
inates, and the electrons are going to be concentrated 
around the center of the sphere anyway. However, the 
greater the radius the less important the confinement, 
and the shape of the XC potential will play a more promi- 
nent role. This trend can be seen in Fig. where we 
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FIG. 4: Exact (solid line) and self-consistent GGA (dashes) 
and EXX+M (dash-dots) densities n (r) and reduced radial 
densities r 2 n(r). All the quantities have been scaled with 
the radius of the confining sphere. The self-consistent LDA 
and MGGA densities are very similar to the GGA ones and 
have not been included in the figure. None of the approx- 
imate functionals is able to reproduce the correct behavior 
of the density in the center of the sphere, and their overall 
performance is very poor for R S> 0. 



compare the DFT densities n{r) as well as the the cor- 
responding reduced radial densities r 2 n(r) with their ex- 
act counterparts. In the atomic- like regime (i? < 5), we 
can observe genuine errors on the DFT densities around 
the center of the sphere. However, this wrong behavior 
will lead to marginal errors on integrated quantities, as 
suggested by the overall good agreement shown by the 
reduced radial densities. At intermediate mean densities 
(R ~ 10), the differences on r 2 n(r) can be already ob- 
served at a first glance. Moreover whereas the system 
shows an incipent localized behavior, which is character- 
ized by a density reaching a local maximum at r =/= 0, 
the EXX+M self-consistent density still has an atomic 
behavior characterized by a maximum at r — 0. Finally, 
in the low density limit, all the approximate functionals 
fail to describe the exact density profile with a minimum 
accuracy. For instance, if R = 25 the error given by (@J 
is 8.6 % using the EXX+M functional and 14.2 % using 
the GGA. 

Then, the full minimizimation of E[n] adds a fur- 
ther source of error due to the inaccuracies of the self- 
consistent density. Nonetheless, these self- consistency- 
induced errors in the total energies are going to be less im- 
portant because there is an overall trend to cancellation, 
see Table [H] For instance, the EXX+M self-consistent 
density is smoother than the exact one, which lowers 
the kinetic and exchange energies while increasing the 
Hartree interaction energy. However, in spite of the dis- 
torted density profile, there is a fortunate cancellation of 



TABLE II: Self-consistent DFT-KS results for several radii 
compared with the exact ones. The last file of each entry 
contains the DFT density error Sn as defined in Eq. 0] 





Exact 


LDA 


GGA 


MGGA 


EXX+M 


R= 1 












-Etot 


11.5910 


11.7338 


11.6376 


11.5540 


11.5770 


T S + Wh 


13.3999 


13.3955 


13.3974 


13.3973 


13.4000 


Ex 


-1.7581 


-1.5193 


-1.6820 


-1.7785 


-1.7582 


Ec 


-0.0507 


-0.1423 


-0.0778 


-0.0648 


-0.0647 


Sn 




0.0070 


0.0049 


0.0056 


0.0006 


R = 5 












Etot 


0.7016 


0.7104 


0.7017 


0.6897 


0.6975 


T s + Wh 


1.0734 


1.0713 


1.0716 


1.0720 


1.0747 


E x 


-0.3335 


-0.2897 


-0.3196 


-0.3401 


-0.3348 


Ec 


-0.0383 


-0.0713 


-0.0504 


-0.0421 


-0.0423 


Sn 




0.0210 


0.0222 


0.0196 


0.0157 


R= 10 












Etot 


0.2381 


0.2371 


0.2346 


0.2305 


0.2362 


T s + Wh 


0.4261 


0.4245 


0.4246 


0.4249 


0.4275 


Ex 


-0.1592 


-0.1395 


-0.1537 


-0.1643 


-0.1608 


Ec 


-0.0288 


-0.0479 


-0.0363 


-0.0301 


-0.0306 


Sn 




0.0383 


0.0448 


0.0371 


0.0401 


R = 25 












Etot 


0.0633 


0.0587 


0.0584 


0.0582 


0.0626 


Ts + Wh 


0.1387 


0.1377 


0.1376 


0.1377 


0.1395 


t^JX 


-u.uoyu 


-u.uooo 


-u.uoy / 


-U.UDoZ 


-U.UDUU 


Ec 


-0.0163 


-0.0256 


-0.0195 


-0.0163 


-0.0170 


Sn 




0.1352 


0.1424 


0.1300 


0.0858 


R = 50 












Etot 


0.0249 


0.0197 


0.0199 


0.0204 


0.0243 


Ts + Wn 


0.0620 


0.0621 


0.0622 


0.0623 


0.0624 


Ex 


-0.0278 


-0.0270 


-0.0318 


-0.0330 


-0.0282 


Ec 


-0.0093 


-0.0154 


-0.0105 


-0.0088 


-0.0099 


Sn 




0.3317 


0.3348 


0.3343 


0.1084 



errors that makes the sum of these three terms practi- 
cally equal to the exact value. Thus, the minimization of 
the energy leads to changes on the density profile favoring 
lower correlation energies, which is done by increasing the 
density around the center of the sphere. As a result, self- 
consistent total energies are just slighty worse than the 
non-self-consistent ones, the relative error AE to t ranging 
from -0.1% for R — 1 to -2.5% for R = 50. A similar con- 
clusion reads for the remaining functionals (LDA, GGA, 
MGGA), although in this case the self-consistent Ts[n] 
is greater than the exact one, whereas the classical inter- 
action energy is reduced after the minimization. Hence, 
self-consistency keeps the fairly good quality of the to- 
tal energies in the atomic limit, where there are minor 
changes in the density profile, and for intermediate and 
low densities the dramatic changes on n(r) induce a few 
percent variation on the total energies that, in any case, 
were already too small. 
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IV. CONCLUSIONS 

In this work we have assessed the quality of some 
of the most popular XC functionals used in ab-initio 
electronic structure calculations in a simple system of 
electrons where several correlation regimes can be easily 
achieved. None of the exchange functionals is completely 
free of self-interaction errors, making imposible to obtain 
accurate energies for highly correlated electrons. Local 
and GGA correlation exhibit a similar drawback, but the 
MGGA correlation shows an excellent performance for all 
the regimes. Nonetheless, neither LDA/GGA nor MGGA 
can describe highly non-local correlation effects that lead 
to a non-trivial behavior of the correlation potential, so 
seriously affecting the quality of the self-consistent den- 
sities. Due to persistent cancellation of trends, the corre- 
sponding changes on the total energy after minimization 
are less important, but this uncontrolled source of error 
might prevent these approximate functionals from having 
full predictive accuracy. 

The overall bad quality of the self-consistent KS po- 
tential also compromises the evaluation of post-self- 
consistency corrections based on more sophisticated 
methods, like Many-Body Perturbation Theory or time- 
dependent DFT |3j], if the wrong effective potential is 
not corrected as well. On the other hand, an accurate 
description of the XC potential is required, for instance, 
when studying neutral excitations in finite systems us- 
ing time-dependent DFT ^^|. As shown recently by 
Delia Sala and Gorling for those systems having an 
HOMO orbital with nodal surfaces, the exact exchange 
potential tends to a constant if going to infinity over a 
set of zero measure directions. This leads to the appear- 
ance of potential barriers that, although could be of mi- 
nor importance when obtaining the self-consistent static 
results, might be essential if a proper description of all 
the unoccupied KS orbitals was required. Here, although 
in a very different context, potential barriers induced by 
the non-locality of the many-body effects in an electron 
system have been observed as well. 



The limitations observed for this family of proto- 
type systems might have relevance for real molecular or 
condensed-matter systems in some cases. Prospective 
tests of the MGGA functional used in this work show an 
excellent performance for a wide variety of typical molec- 
ular and solid state systems which, moreover, seems to 
be kept if partial self-consistency is achieved [8|. How- 
ever, there is no guarantee that the energy minimization 
procedure may lead to small, but relevant changes on the 
local electron density in, for instance, the bonding region 
between two species, so giving a wrong account of the 
nature of such a chemical bond. Finally, for this model 
system we have seen that there are not substantial differ- 
ences between the fully self-consistent GGA and MGGA 
densities, although in this case the MGGA functionals 
take a simple GGA-like form. In spite of this simplifi- 
cation, the MGGA-XC potential amplifies the spurious 
oscillations already appearing under the GGA prescrip- 
tion. As a conclusion, the overall capability of MGGA 
and envisaged improvements thereof to yield accurate XC 
energies can hardly be questioned (specially for the cor- 
relation part), but further studies of the actual perfor- 
mance of the corresponding KS potentials are required. 
The latter point is relevant for those situations in which 
LDA/GGA are not able to reproduce the electron density 
with the requiered accuracy, and for other DFT-based ap- 
plications needing an overall good description of «s(r). 
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